LCAO model for 3D Fermi surface of 
high-T c cuprate Tl 2 Ba 2 Cu0 6+( 5 



ON 
O 

o 

Or 
<D 
00 

in 



o 
o 

Oh- 



o 

> 

o 

On 

o 
in 
o 

-I— > 
c3 



i 

c 

o 
o 



X 



M. Stoe^E and T.M. Mishono\§ 

Department of Theoretical Physics, Faculty of Physics, 
University of Sofia St. Kliment Ohridski, 
5 J. Bourchier Boulevard, BG-1164 Sofia, Bulgaria 
(Dated: September 15, 2009) 

A simple analytical formula for three-dimensional Fermi surface (3D FS) of T^I^CuOe+s is 
derived in the framework of LCAO approximation spanned over Cu 4s, Cu 3d x 2_ y 2 , O 2p x and 2p y 
states. This analytical result can be used for fitting of experimental data for 3D FS such as polar 
angle magnetoresistance oscillation. The model takes into account effective copper-copper hopping 
amplitude t ss between Cu 4s orbitals from neighbouring Cu02 layers. The acceptable correspondence 
with the experimental data gives a hint that the t as amplitude dominates in formation of coherent 
3D FS, and other oxygen-oxygen and copper-oxygen amplitudes are rather negligible. For absolute 
determination of the hopping parameters a simple electronic experiment with a field effect transistor 
type microstructure is suggested. The thin superconductor layer is the source-drain channel of the 
layered structure where an AC current is applied. 

PACS numbers: 71.15.Ap, 73.43.Qt, 74.72.Jt, 74.78.Fk 



I. INTRODUCTION AND NOTATIONS 

For a long time magnetic oscillations are a stan- 
dard method for determination of Fermi surface. For a 
comprehensive introduction see for example the mono- 
graph by Shoenberg*. A recent observation of three- 
dimensional Fermi surface (3D FS) in T^E^CuOg+s 
(T1:220T)2 unambiguously has shown that charge carri- 
ers in this material are ordinary Fermi particles. This 
observation has an important significance for the physics 
of high-T c cuprates in general. There is almost a consen- 
sus that pairing mechanism is common for all cuprates 
wherever it is hidden. That is why the observation of 3D 
FS leads to the conclusionthat this mechanism should be 
able to work even for Fermi quasiparticles in BCS sce- 
nario. An important first step in this scenario is the 
analysis of the electron spectrum of a metal in a selfcon- 
sistcnt approximation of independent electrons. For the 
cuprates, as for many other ion materials, the method of 
linear combination of atomic orbital (LCAO) gives ade- 
quate description of the band structure. Tight binding 
band structure is also a relevant starting point for theo- 
retical analysis of many phenomena related to quasiparti- 
cle interaction. For example Abrikosov^ used tight bind- 
ing cr-model to build the spin density theory of metal- 
insulator transition in cuprates. A lot of phenomena 
in cuprate physics especially for overdoped cuprates can 
be understood in the framework of orthodox fermiology. 
The purpose of the present work is to derive analyti- 
cal LCAO formula for 3D FS and the band structure of 
Tl:2201 which can be used for further analysis of the ex- 
perimental data; for local density approximation (LDA) 
calculations of band structure of this material see, for 
example, reference^. These analytical results can be 
used for fitting the angle magnetoresistance oscillation 
(AMRO) data and at the same time they represent a re- 
alistic noninteracting part of the lattice Hamiltonian for 



further consideration of pairing in cuprates. 

Hilbert space in LCAO approximation is spanned over 
the relevant Cu 3d x 2_ y 2, Cu 4s, O 2p x and O 2p y or- 
bitals. The generic 4-band cr-model for Cu02 plane was 
suggested by Labbe and Bok— . Later on detailed ah ini- 
tio calculations of band structure of layered cuprates by 
Andersen et al~ confirmed that this generic model ad- 
equately interpolates the LDA band structure and the 
influence of 7r-orbitals for the conduction band is neg- 
ligible. In TL2201 the neighbouring Cu02 planes are 
shifted in a half period. In an elementary cell indexed by 
three integer numbers R n = aini + a 2 n 2 + 2*3113, where 
ai = a e x , a 2 = a e y and a 3 = 6 e z + \a (e x + e y ), the 
space vector of copper ions is Rc u = 0, and for oxygen 
ions we have Ro a = \aoe x and Ro b = \aoe y ; e x , e y and 
e 2 are the unit vectors. For the introduced notations the 
LCAO wave function reads as 



^LCAo(r) = ErjAiV'cuadO R n RCu) 
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where * n = (D n , S n , X n ,Y n ) is the tight-binding wave 
function in lattice representation^. In second quantisa- 
tion approach the complex amplitudes D n , S n , X n and 
Y" n become Fermi annihilation operators. 

The LCAO wave function can be expressed by Fermi 
operators in momentum space 
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where Af is the number of unit cells supposing peri- 
odic boundary condition, k = (j> x /ao,p y /ao,p z /bo) is 
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the electron quasimomentum and p is the dimension- 
less quasimomentum in first Brillouin zone (p x ,p y ,p z ) G 
(— 7T, it). This equation describes the Fourier transfor- 
mation between the coordinate representation f n = 
(D n , S n , X n , Y n ) and the momentum representations 
ip p = (D p , S p , X p , Y p ) of the tight-binding wave function. 

In order to derive the electron band Hamiltonian in 
momentum representation, we have to analyze the LCAO 
Schrodinger equation etpi = e^i - Y2j(^ z )^ij 1 Pj f° r a pl ane 
wave ^> n oc e lkRn , where ipi are amplitudes multiply- 
ing atomic orbitals, e is the electron energy, ej are sin- 
gle site energies, the sum is spanned to the nearest and 
next nearest neighbouring atoms at i-orbital and hop- 
ping amplitudes have +sign for bonding, and -sign for 
antibonding orbitals with different signs of atomic wave 
functions^. For a detailed analysis, notations and refer- 
ences seem. Here we will give a brief description. Let 
e(p) be the electron band dispersion; ed is the single-site 
energies for Cu 3d level, t pd is the hopping amplitude 
between the O 2p states and Cu 3d. Considering only 
nearest-neighbour hoppings and omitting the common 
gjk-Rn mu itiplier for the Cu 3d x 2_ y 2 amplitude in R n 
elementary cell we obtain 

e(p)Dp = e d Dp-it pd (e tE ?X p -e tE ?e- l P*X p -e^Yp+ 
e^e-^Yp) = e d D p + 2t pd (sm £fX p - sin p -fY p ) = 

(3) 

Analogous consideration can be applied for Cu 4s ampli- 
tude: e s is the single site energy for Cu 4s state, t ap is the 
hopping amplitude between O 2p and Cu 4s; t sp is bigger 
than t p d because Cu 3d x 2_ y 2 orbitals are much more lo- 
calized near to copper nucleus. Due to the symmetry of 
the Cu 4s wave function, now the relative signs between 
X p and Y p are equal 

6(p)*Sp — £s*Sp tspSxXp -\- t S pSyY p 

-\-it ss (e lS ?~ e 1 ^ e lp * + e 1 ^ e 1 ^ e~ Wz + e~ l2 r e lVz 

■ p x ■ Py ■ ■ p x ■ Py ■ ■ p x ■ Py 

_l_ e -»-2- e »— e -ipz + e ~ l — e -i— e Wz _|_ e - *"2~e -i ~2" e~ lPz 
+ e ^ e -^ e ip. +e iE Te~ lE Te-* p *)Sp = 

where z(p) = c x c y c z . The influence of three dimension- 
ality (3D) is taken into account only here by the effective 
Cu 4s-Cu 4s transfer integral t ss between copper ions 
from different Cu02 planes. Following!* we have used the 
standard notations: 

s x = 2sin(p x /2), s y = 2sin(p y /2), 
c x =2cos(p x /2),c y = 2cos(p y /2), (5) 
x = sin 2 (p x /2), y = sin 2 (p y /2), 



adding c z = 2cosp 2 . In first Brillouin zone for p x ,p y £ 
(— 7r,7r) the variables c x ,c y > 0, however in the whole 
momentum space in the interval (0, 2ir), for example, we 
have to redefine c x — 2\ cos(p x /2)\ and c y — 2| cos(p x /2)\. 

Analogously considering electron hopping to O 2p.,; or- 
bital and dividing by —ie 1- ^ we have 

e (p)A p = e p Ap + t pd s x D p + t sp s x S p + 

+t pp (— i\e 2 e 2e * a e a +e » e a — e a e 2 )ip 

— £pXp -\- tpdS x Dp -\- t S pS x Sp tppS x SyYp, 

(6) 

where e p is O 2p single-site energy, t pp is the hopping 
amplitude between adjacent O 2\i x and O 2p y orbitals. 
Due to the crystal symmetry the equation for Y p can be 
obtained by exchange between x and y, and X and Y; 
only the relative sign between O 2p and Cu 3d orbitals 
has to be changed. 

Finally the LCAO Schrodinger equation in momentum 
representation reads as 

(Hp - e(p)l)^ = 0, (7) 

where 
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(8) 

and ip p = (Dp, S p , Xp,Y p ) is the wave function in mo- 
mentum space. 

After some algebra the secular equation takes the form 

D(p) =det(H p -e(p)l)=^xy + 6(x + y)+C 

(9) 

+ z[Kxy + C(x + y) + M] = 
with energy-dependent coefficients 
A(e) = 16(4i 2 d i s 2 p + 2i 2 p i pp e d - 2i 2 d i pp e s - t 2 p p e d s s ), 

B(e) = -Ae p (t 2 sp e d + t 2 pd e s ), C(e) = e d e 2 p e s , 

JC(e) = ~t ss t pp (e d tpp + 2i 2 d ), £(e) = -e p £ ss i 2 d , 

M(e) = e d elt ss , 

(10) 

where e s = e(p)-e B , £ P = e(p)-e p , e d = e(p)-e d , are the 
energies taken into account from single site atomic levels. 
Due to the small numerical value of t ss the modulation 
of FS in p z direction is also small. 



II. ANALYSIS OF THE INFLUENCE OF THE 
INTERLAYER HOPPING AMPLITUDE t« 

Comparing our secular determinant with the purely 
2D case t ss = 0, in ^ one can see that the influence of 



3 



interlayer hopping is formally reduced to a momentum 
dependent single site energy shift for the Cu 4s level 



t ss z(p). 



(11) 



This is a diagonal matrix element whose influence is just 
zero on the p y = tt and p x — tt lines in 2D Brillouin zone, 
shown in figure 1, which we will discuss later. On the 
diagonal, where p x = p y and p y — 2tt — p x , the influence 
of this term is also negligible. The 2D Hamiltonian (for 
t S s = 0) has such eigenfunctionsii: 
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(12) 

After the calculation of ijj p 
(p) we have to make the normalization 



= i z - 

c(2D) 



2 tst PP- 



where 
for e : 

ip.-p := ifjp/ Dp + Sp + X 2 + Y 2 . In first perturbational 

approximation the influence of the diagonal term (111[) 
gives an addition to the band energy e(p) 



W(p) = -8t ss cos (f ) cos (&l ) cos( Pz )S 2 p 
s(p)=s^ D Hp) + W(p), 



(13) 



i.e we present the band energy as a sum of a 2D band 
energy and a correction by taking into account the hy- 
bridization in c-direction W(p). Since S p cx (x — y) this 
correction reduces to zero on the diagonals of Brillouin 
zone and in such a way the influence of t sa vanishes at 8 
points of 2D Fermi contour at every 45 degrees at hori- 
zontal, vertical and diagonal lines crossing (it, tt) point. 
For the overdoped Tl:2201 the hole pocket centered at 
(tt, tt) point takes 62% of the Brillouin zone, see figurel. 

In initial approximation t ss — we have a 2D secular 
equation 

det(H p - e ( 2D > (p)l) = .4xy + B(x + y) + C, (14) 

which gives explicit expressions for the lower and upper 
arch of constant energy contour (CEC)^ 



Py OV,) (p x ) — 2 arcsin J— 



Bx+C 
Ax+B 1 



(15) 



Pv V \Px) = 27T -p y (px). 



For the band velocity at energy equal to Fermi energy 
e^(p)=E F 

_ 8fS 2D \p) 



v F = \lvl + vl, (16) 
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we also have simple explicit expressions 10 

_ 1 (Ay+B) sin( P:c ) 
Ux ~ 2 A'xy+B'(x+y)+C" 

_ 1 (Ax+B)sm(p y ) 
U V — 2 A'xy+B'(x+y)+C ' 



(17) 



dp 



Hp) 



_ [(Ay+B) 2 x(l~x) + {Ax+Bfy{l-y)] in _ 
~ \A'xy+B'{x+y)+C'\ ' 



(18) 



the velocity in m/s is actually ao~v/h. The coefficients in 
the denominator A' , B' and C are energy derivatives of 
the polynomials (JTDJ) , 

^'(e) = 16 [2i s 2 p i pp - 2t 2 d t pp - t 2 pp (e d + e 8 )] = dA/de, 
B'(e) = -4(t 2 e d + t 2 e s ) - Ae p (t 2 + t 2 ) = dB/de, 



C'(e) = e s e 2 + e d £p + 2e d e s £ P = dC/de. 

(19) 

Our problem is to take into account the influence of the 
perturbation (fT3"|) to the Fermi contour (fT5]) . Under the 
influence of the perturbation W^p) every point p of this 
CEC is shifted in a perpendicular to CEC direction with 
momentum shift 



Ap = -W(p)£, «* = W' 
e(p) = £ (2D) (p)+V'Ap. 



(20) 



The coordinates of the perturbed CEC are p x + Ap x and 
p y + Ap y ; in such a way we approximatively built the 
3D Fermi surface. The formulae above represent a self- 
explainable derivation: (1) the band energy is approxi- 
mately presented by gradient expansion (2) the ratio of 
energy difference and momentum difference is equal to 
Fermi velocity (3) the shift of Fermi contour in momen- 
tum space is parallel to the Fermi velocity. Projections 
of this Fermi surfaces in 2D Brillouin zone are depicted 
in figure 1. The LCAO approximation gives a similar 
shape of the Fermi contours as the experimental Fermi 
surface^. We have taken the set of parameters from 6 : 
e s = 6.5 eV, e d = eV, e p = -0.9 cV, < pd = 1.6 eV, 
i pp = eV, t sp = 2.3 eV. The Fermi level E F = 1.89 eV 
is determined to give / = 62% hole filling 
2 of the 2D Brillouin zone 



(2^) 2 



(Vx ~ Py(Px] E F )) dpx = f, 



(21) 



where pd is the solution of the equation < p y (pd] E F ) 
Pd < tt- 



Xd = sin 



(y) =j(-B+Vb 2 ~~Ac). (22) 



The Fermi contour pases trough the points D — (pd, Pd) 
and C — (p c , tt), where 



x n = sin 



• 2 (Pc\ 

m It) 



p c \ B + C 
' A + B' 



(23) 



or p c = p v (tt; E F ). Finally interlayer hopping amplitude 
t ss = 140 meV is determined by the comparison with the 
experiment of the modulation of FS 2 . 
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FIG. 1: Projections of 3D FS in 2D Brillouin zone obtained 
by equation (|20|) . In agreement with AMRO data£ the small 
modulation due to Cu 4s-Cu 4s tunnelling amplitude t ss van- 
ishes at 8 highly symmetric lines passing through (n, n) point. 
The hole pocket centered at (ir, ir) point has approximately 
62% of the surface of the 2D Brillouin zone. The warping is 
not to scale, it is exaggerated aftei^ in order to emphasise its 
angular dependence which can be explained in the framework 
of the LCAO method. 



III. COMPARISON WITH ARPES DATA 

The galvanomagnetic phenomena such as AMRO are 
sensitive mainly to the total area of the sections of 
the Fermi surface. The angle resolved photoemission 
spectra (ARPES) however are sensitive to the shape 
of the quasi 2D Fermi surface^. In order to make a 
compromise conserving the area of the Fermi surface, 
actually on the cross-section at p z = 0, we can try 
to fit its shape. The Fermi contour extracted from 
ARPES data for TL2201— we can use the diagonal 
point D = (0.35767r/ao,0.35767r/ao) and another point 
C = (0.1256 7r/ao,7r/ao) as reference points. We can start 
from the set of LCAO parameters given in reference^ and 
changing only the Fermi level Ep and Cu 4s level e s we 
can pass the Fermi contour through the reference points 
C and T> as it is done in figure 2. The Fermi contour re- 
produces the shape from reference^, its area is in agree- 
ment with the AMRO dataS, and even the Ep and e s 
are not very different from the fit of LDA calculations by 
Andersen et al&. 

The fit of the absolute value of the energy scale however 
requires a big compromise. As it is well-known the LDA 
often gives overbinding of order 2 or even 3. Correcting 
overbinding in local-density-approximation calculations 
we can insert an energy renormalization scale so that the 
shape of the Fermi surface to be exactly conserved and 
only the energy width along some well investigated cut to 
coincide with the experiment. We used the cut-III from 
reference^ to fix the energy scale. This cut is given as a 
short segment in figure 2. In figure 3 this cut, energy ver- 
sus quasimomentum, is given as leftmost segment. The 
circles trough which the dispersion line passes are refer- 
ence points chosen from the experimental data^i These 




FIG. 2: A two dimensional section of Fermi surface for 
p z — 0. The theoretical Fermi contour is passing through 
the reference points in fc-space: C = (0.1256 n/ao, 7r/ao) and 
T> — (0.3576 7r/oo, 0.3576 7r/ao) marked with small (red) cir- 
cles. The short segment close to the saddle point (7r/ao, 0/ao) 
is the cut-III from Ref. 13 analyzed further at figure 3. 
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FIG. 3: Energy (in arbitrary units) as function of momentum 
along the cut-ffli 3 " presented in quasimomentum space in fig- 
ure 2. The two (red in the color version) circles are reference 
points determining the vertical shift of the dispersion curve 
and the renormalization of all LDA energy parameters. For 
comparison are presented standard theoretical cuts e(p) for 
the triangle: (0,0)-(7r,0)-(7r,7r)-(0,0). 



reference points determine the energy scale. 

On the right of this reference segment is presented the 
standard energy dispersion along highly symmetric direc- 
tions (0,0)-(7r,0)-(7r,7r)-(0,0) in momentum space. One 
of the reference points was shifted in vertical direction 
in order to use only the energy width of the cut-III for 
energy determination. We believe that the so determined 
band structure has predictive capabilities in a sense that 
we can use analytical formulae for the band structure to- 
gether with the fitted values of the LCAO parameters 
in order to predict the results of new ARPES measure- 
ments. In short, the mission of the theoretical physics 
is to predict the results of experiments not done by no- 
body before. As a further perspective we believe that 
single site energies have to be taken from the experimen- 
tal data processing of spectroscopic data for interband 
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FIG. 4: A field effect transistor(FET) is schematically illus- 
trated. The current I(t), applied between the source (S) and 
the drain (D) has a frequency u. Running through the tran- 
sistor the electrons create voltage Usg with double frequency 
2uj between the source (S) and the gate (G). The source-drain 
voltage Usd is measured on the triple frequency 3uj. 



transitions. Of course such an interpretation suppose ac- 
tive use of the generic 4-band model taking into account 
the accessories for every cuprate. Every additional exper- 
iment for absolute determination of the energy scale can 
be an indispensable tool for the final determination of the 
energy scale of LCAO parameters. In the next section we 
propose a simple electronic measurement which should 
be done with a layered metal-insulator-superconductor 
structure with the same superconductor. 



IV. DETERMINATION OF LOGARITHMIC 
DERIVATIVE OF DENSITY OF STATES BY 
ELECTRONIC MEASUREMENTS 

The analysis of Fermi surface by AMRO and ARPES 
as a rule leads only to the determination of relative val- 
ues of the parameters of the LCAO Hamiltonian. In or- 
der to extract the absolute values we need to fit band 
widths, cyclotron frequencies, density of states v{Ep) at 
Fermi energy Ep or other quantities having dimension of 
energy. We consider as very important the comparison 
of different methods for investigations of band structure. 
The AMRO gives the exact value of the area of sections 
of Fermi surface, the ARPES is more sensitive for the 
shape and will be nice if this knowledge is completed by 
some other method giving the exact value of some energy 
parameter. 

The purpose of this section is to suggest a simple elec- 
tronic experiment, determining the logarithmic deriva- 
tive of the density of states (DOS) by electronic mea- 
surements using a thin film of the investigated material 
Tl:2201. The proposed experiment requires the prepa- 
ration of field-effect transistor (FET) type microstruc- 
ture and requires standard electronic measurements. The 
FET controls the current between two points but does so 
differently than the bipolar transistor. The FET relies 
on an electric field to control the shape and hence the 
conductivity of a "channel" in a semiconductor material. 
The shape of the conducting channel in a FET is altered 
when a potential difference is applied to the gate termi- 
nal (potential relative to either source or drain) . It causes 



the electrons flow to change it's width and thus controls 
the voltage between the source and the drain. If the neg- 
ative voltage applied to the gate is high enough, it can 
remove all the electrons from the gate and thus close the 
conductive channel in which the electrons flow. Thus the 
FET is blocked. 

The system, considered in this work is in hydrody- 
namic regime, which means low frequency regime where 
the temperature of the superconducting film adiabati- 
cally follows the dissipated Ohmic power. All working 
frequencies of the lock-ins, say up to 100 kHz, are actu- 
ally low enough. The investigations of superconducting 
bolometers show that only in the MHz range it is nec- 
essary to take into account the heat capacity of the su- 
perconducting film. As an example there is a publication 
corresponding to this topio^, as well as the references 
therein. In this work we propose an experiment with 
a FET, for which we need to measure the second har- 
monic of the source-gate voltage and the third harmonic 
of the source-drain voltage. Other higher harmonics will 
be present in the measurements (e.g. from the leads), but 
in principle they can be also used for determination of the 
density of states. An analogous experimental research 
has been already performed for investigation of thermal 
interface resistance^ The suggested experiment can be 
done using practically the same experimental setup, only 
the gate electrodes should be added to the protected by 
insulator layer superconducting films. 

The purpose of this section is to suggest a simple elec- 
tronic experiment, determining the logarithmic deriva- 
tive of the density of states by electronic measurements 
using a thin film of the material TL2201. The thickness 
of the samples should be typical for the investigation of 
high-Tc films, say 50-200 nm. Such films demonstrate 
already the properties of the bulk phase. The numerical 
value of this parameter 



<M«0 
de 



(24) 



will ensure the absolute determination of the hopping 
integrals. 

We propose a field effect transistor (FET) from Tl:2201 
Fig. [4] to be investigated electronically with a lock-in at 
second and third harmonics. Imagine a strip of Tl:2201 
and between the ends of the strip, between the source (S) 
and the drain (D) an AC current is applied 



hr>(t) = I cos(u>t). 



(25) 



For low enough frequencies the ohmic power P in- 
creases the temperature of the film T above the ambient 
temperature To 



P = Rll D =a(T-T ) 1 



(26) 



where the constant a determines the boundary thermo- 
resistance between the Tl:2201 film and the substrate, 
and R(T) is the temperature dependent source-drain 
(SD) resistance. We suppose that for a thin film the 
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temperature is almost homogeneous across the thickness 
of the film. In such a way we obtain for the temperature 
oscillations 

( T _ To ) = M^ = M cos 2 (6L ,i ) . (27 ) 
a a 

As the resistance is weakly temperature dependent 

dR(T) 



where the logarithmic derivative of the density of states 
v(t) taken for the Fermi energy Ep has dimension of 
inverse energy, the work function W has dimension of 
voltage, T is the temperature in Kelvins and ks is 
the Boltzmann constant. For an introduction see the 
standard textbooks on statistical physics and physics of 
metals . 25 ' 26 Substituting here the temperature variations 
from Eq. gives 



R(T) = R Q + (T-T )Ro ', Rq '(T ) 



dT 



(28) 
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71-2 k l V 'rra 

' 6e v 



1 



Rplq 
aT 



cos(2wt) 



O(/ 4 ), (36) 



A substitution here of the temperature oscillations from 
Eq. I|2?p gives a small time variations of the resistance 



R(t) = R ( 1 + ^/2 cos 2 (wi) 
a 



Now we can calculate the source-drain voltage as 



U SD (t) = R(t)I SD (t). 



(29) 



(30) 



where the O-function again marks that the terms having 
Iq are negligible. 

The oscillations of the temperature creates AC oscilla- 
tions of the source-gate (SG) voltage. We suppose that a 
lock- in with a preamplifier, having high enough internal 
resistance is switched between the source and the gate. 
In these conditions the second harmonics of the work- 
function and of the SG voltage are equal 



Substituting here the SD current from Eq. (|25|) and the 
SD resistance from Eq. (|2"9")) gives us for the SD voltage 
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7T A kf 



SG 



"^0 aT ■ 



U SD {t) = U^ f) cos(cjt) 



u, 



f D f) cos(3urf). (31) tfsoW = cos(2 w t) + Lf G /} cos(4^) + 



(37) 



The coefficient in front of the first harmonic f/gJ/' ) « RqIq 
is determined by the SD resistance Rq at low currents Iq, 
while for the third harmonic signal using the elementary 
formula cos 3 (ut) = (3 cos (uit) + cos (3wt))/4 we obtain 



U. 



(3/) 
SD 



U. 



(1/) 
SD 



4o 



IqR'q. 



(32) 



From this formula we can express the boundary thermo- 
resistance by electronic measurements 



U, 



(i/) 



SD 



4CA 



(3f) 



(33) 



SD 



The realization of the method requires fitting of R(T) 
and numerical differentiation at working temperature To; 
the linear regression is probably the simplest method if 
we need to know only one point. 

At known a we can express the time oscillations of the 
temperature substituting in Eq. (12711 



Substituting a from Eq. (|33|) we have 



(2/) _ 



SG 



Ge 



SD 

Jo 



R' 



(38) 



From this equation we can finally express the pursued 
logarithmic derivative of the density of states 



dlnz/(e) 



de 



E p 



3e I Q U^J ] dR 

2^fc|^^WTdr- 



(39) 



In such a way the logarithmic derivative of the density 
of states can be determined by fully electronic measure- 
ments with a FET. This important energy parameter can 
be used for absolute determination of the hopping inte- 
grals in the generic LCAO model. The realization of the 
experiment can be considered as continuation of already 
published detailed theoretical and experimental investi- 
gations and having a set of complementary researches we 
can reliably determine the LCAO parameters. 



T = To + ^ (1 + cos(2^)) w To ( 1 + cos ( 2wt ) 
2a \ 2aT 

(34) 

In this approximation terms containing Iq are neglected 
and also we consider that the shift of the average tem- 
perature of the film is small. 

The variations of the temperature lead to a variation of 
the work function of the film according to the well-known 
formula from the physics of metals 



6e v 



(35) 



V. DISCUSSION AND CONCLUSIONS 

The suggested LCAO model for FS of Tl:2201 
describes the important qualitative properties of p z 
modulation: vanishing of this modulation at 8 sym- 
metric points of the 2D Fermi contours. This is an 
important hint that t ss amplitude dominates in the 
formation of coherent 3D Fermi surface and other 
interlayer tunnelling amplitudes are irrelevant. This 
qualitative conclusion for the importance of Cu 4s c-axis 
tunnelling matrix elements is in agreement with the 
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long-predicted analysis by Andersen et al£. We wish to 
add a few words concerning the pairing mechanism in 
Tl:2201 and cuprates in general. Wherever it is hidden 
its influence has to be reduced to effective momentum 
dependent scattering amplitude for electron pairs from 
conduction 3d band with opposite momcntums; band 
structure created by 3d-2p-4s hybridization for which 
the revealing of FS is an indispensable first step. 
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VI. APPENDIX A: CALCULATION OF 
EIGENVECTOR 



In the equation for eigenvectors (Hp 2D ^ — e(p)l)^i p , 

with Hp 2D ^ from f8j we search for a solution in the form 
ip p = (1, S p , X p , Y p ). In such a way we obtain the system 



-e s 




tspSy \ 






-Bp 


— tpp^xSy 1 t 






tpp^xSy 




V ^ 




The solution 5 = A s /A, X = A x /A, Y = A Y /A is 
presented by the determinants 



det 



-t 



det t 



A Y = det 



e s 

sp^x 



L sp c> x 
pp-^a^y 



o 

"^pd'Sa; 



^spSx 



L sp°x 



^spSy 
~tpp$x$y 
-En 



~tppS x Sy 




(42) 



(43) 



(44) 



det 



^sp^ic 
tspSy 



L sp°x 

-£ P 

tppSx Sy 



^spSy 
'tppSgSy 
-En 



(41) 



Multiplication by A gives the eigen vector (D, 5, X, V) = 
(A, As, Ax, Ay) presented in equation (fl"2ll . 



